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Abstract. Strongly interacting, dilute Fermi gases exhibit a scale-invariant, universal 
thermodynamic behaviour. This is notoriously difficult to understand theoretically 
because of the absence of a small interaction parameter. Here we present a systematic 
comparison of theoretical predictions from different quantum many-body theories with 
recent experimental data of Nascimbne et. al. [Nature 463, 1057 (2010)]. Our 
comparisons have no adjustable parameters, either theoretically or experimentally. All 
the model approximations seem to be fluctuating around and not converging towards 
the experimental data. It turns out that a simple Gaussian pair fluctuation theory gives 
the best quantitative agreement, except at the critical superfluid transition region. In 
the normal state, we also calculate the equation of state by using a quantum cluster 
expansion theory and explore in detail its applicability to low temperatures. Using the 
accurate experimental result for the thermodynamic function >5'(T), we determine the 
temperature T of a trapped Fermi gas at unitarity as a function of a non-interacting 
temperature Ti which can be obtained by an adiabatic sweep to the free gas limit. 
By analyzing the recent experimental data, we find a characteristic temperature 
(T/Tf)o = 0.19 ± 0.02 or {T,/Tf)o = 0.16 ± 0.02 in a harmonic trap, below which 
there are deviations from normal Fermi-liquid-like behavior that may be attributed to 
pairing effects. Here Tp is the Fermi temperature for a trapped ideal, non-interacting 
Fermi gas. Our thorough comparison may shed light on further theoretical development 
of strongly interacting fermions. 
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1. Introduction 

The recent discovery of broad Feslibach resonances in two-component atomic Fermi gases 
has opened a new era in the study of strongly interacting fermions [B [21 |3l HI O |6l [7] . 
By tuning an external magnetic field across the Feshbach resonance, the interatomic 
attractions can be changed precisely from weak to infinitely strong, leading to the 
observation of a crossover from a Bardeen-Cooper-Schrieffer (BCS) superfluid to a Bose- 
Einstein condensation (BEC). At resonance, the s-wave scattering length diverges 
(a<, = ±oo) and the two-body scattering amplitude reaches the maximum value allowed 
by quantum mechanics due to unitarity. Many unique properties are anticipated in this 
strongly interacting limit, including a high superfluid transition temperature and an 
exotic normal state with a pseudogap. 

Most interesting of all is fermionic universality. This means that all strongly 
interacting, dilute Fermi gases behave identically, regardless of the details of the 
interaction. Their properties depend only on temperature, together with a scaling factor 
equal to the average particle separation [SI [9l [TOl [11] . This limit promises to bring a new 
rigour and simplicity to the understanding of strongly correlated Fermi gases. Because 
of universality, it should be feasible to understand other strongly interacting Fermi 
superfiuids from experiments in the highly controlled environment of an atomic physics 
laboratory. Possible examples include neutron stars and high-Tc superconductors. 

Intense experimental investigations have been carried out to understand fermionic 
universality, in particular, its implication to the thermodynamic properties [121 [131 [HI 
[T5l [T6l [T71 [T8l [T9l [20] . Pioneering observations were carried out at Duke University 
[T2l [T3l [HI [15] , which made the first attempt to reach the unitarity limit in ®Li gas in 
2002 [12]. The stability of atomic Fermi gases with strongly attractive interactions was 
observed. The ground state energy was found to be reduced significantly compared to 
its ideal, non-interacting limit. The reduction factor /3 (or ^ = 1 + /3) has now been 
determined accurately to within a few percent: /3 ~ —0.59 ± 0.01, after substantial 
experimental effort |20] . 

In this paper we show that recent highly accurate measurements [19] on strongly 
interacting ^Li give the most stringent test to date of fermionic strongly coupling 
many-body theories. In fact, these experiments determine the whole set of universal 
thermodynamic functions for a trapped Fermi gas at unitarity. As well as measuring bulk 
thermodynamic properties, the data can be used to determine the energy and entropy, 
E{T) and S(T), of a trapped gas. Due to the use of larger samples, the accuracy is even 
better than that achieved at Duke. With unprecedented precision, these new universal 
functions therefore provide an unbiased test of theoretical predictions. This comparison, 
without any fitting parameters, indicates that while a BCS-type mean field theory is 
certainly incorrect, an extension using a simple Gaussian pair fluctuation theory provides 
overall the best agreement with experiment, except at the critical superfluid transition 
region. 

To give some background to these developments, the first energy and heat capacity 
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measurements as a function of temperature were performed by Kinast et al. |13j . 
However, due to the lack of reliable thermometry in the strongly interacting regime, an 
empirical thermometry was used. Conversion of measured results to real temperature 
required a particular strong-coupling theory, and was therefore mo del- dependent. This 
difficulty of model dependence was overcome at the end of year 2006 [HI [15], by means 
of direct measurements of entropy instead of temperature. In this way, both energy E 
and entropy S were determined without invoking any specific theoretical model. These 
pioneering and very important model-independent measurements had an accuracy at 
the level of only a few percent. 

At the same period, the potential energy of a strongly interacting ^''K gas was 
also measured at JILA fTTJ. The temperature was characterized in terms of the non- 
interacting temperature of an adiabatically equivalent ideal Fermi gas. This is therefore 
equivalent to an entropy measurement. These sets of experimental data, together with 
results of another ^Li experiment at Rice, were analyzed by the present authors [11]. The 
result was that all the thermodynamic data lay on a single universal curve. This gave 
the first, very strong evidence for the universal thermodynamics of a strongly interacting 
Fermi gas [TT] . 

In parallel with these ground-breaking experiments, there have been numerous 
theoretical studies of the thermodynamics of a strongly interacting Fermi gas. In the 
absence of exact solutions, the methods used were either strong-couphng perturbation 
theories [2ll[22l[23l[2a[2a[26l|271[28l[29l|3a 

ab-imtio quantum Monte Carlo (QMC) methods [lIlll2lll3lll3ll5lll6lll71llHlll9l[50]. 

However, a deeper understanding is made difficult by the absence of a controllable small 
interaction parameter [51]. The use of standard perturbation theories thus requires 
infinite order expansions. These typically require truncations of sets of diagrams 
which cannot be fully justified a priori. Numerically exact QMC calculations are very 
helpful and can provide unbiased benchmarks, provided that there is an appropriate 
extrapolation of the lattice results to large lattice size or, in the diagrammatic Monte 
Carlo case, to zero range potentials. 

The first theoretical explanation of the heat capacity at unitarity was given by 
Chen et al. [13], using a pseudogap theory. In this study an empirical temperature was 
converted to an approximate real temperature. The present authors subsequently gave 
a theoretical prediction for the both the homogeneous and trapped equation of state at 
unitarity [30]. This used a Gaussian pair fiuctuation (GPF) theory below threshold, thus 
extending an approach proposed initially by Nozires and Schmitt-Rink (NSR) [21], [29] 
for the above threshold case. We showed that the conversion of empirical temperature 
to actual temperature is strongly model-dependent. 

Thus, in principle one cannot obtain accurate information about the real 
temperature from these empirical temperature measurements, without a reliable strong- 
coupling theory. The model independent measurement of energy as a function of entropy, 
E{S), by Luo et al. [TU [15] was therefore a crucial experimental advance. This 
provided the first data that could be used to quantitatively compare different strong- 
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coupling theories without any free parameters. One such comparison was performed by 
the present authors [38], by using different perturbation theories and available QMC 
results [lEj. Even so, it was still impossible to determine the dependence of the energy 
on temperature E(T) and of the entropy on temperature S(T), due to difficulties in 
determining the absolute temperature T. Moreover, the measurements and comparisons 
were restricted to the case of a trapped Fermi gas. 

Most recently, a general method was developed by Nascimbne et al. at ENS to 
measure the bulk equation of state of a homogeneous Fermi gas of lithium-6 atoms |19j . 
following a theoretical proposal by Ho and Zhou [52]. The local pressure P{fi{z),T) or 
the local thermodynamic potential Q{fi{z),T) = —P{fi{z),T)6V {6V is the volume of 
a cell at the position z) of the trapped gas was directly probed using in situ images 
of the doubly-integrated density profiles along the long 2-axis. The temperature was 
then determined by using a new thermometry approach employing a ^Li impurity. The 
chemical potential could also be determined using the local density approximation, 
with fi{z) = fio — Vtrap{z) and the central chemical potential /^o being determined 
appropriately. By introducing a universal function [53] 

experimentalists were able to determine h{() with very low noise. Here, ( = 
exp {—fi/kBT), T) is the interacting thermodynamic potential and T) is 

the thermodynamic potential of an ideal two-component Fermi gas. This precise 
measurement allows a direct comparison with many-body theories developed for a 
uniform Fermi gas [T9] . 

Our main results may be summarized as follows. First, though the theoretical 
predictions from all the model approximations seem to fluctuate around the experimental 
data, it turns out that the simplest Gaussian pair fluctuation theory gives the best 
description of the observed thermodynamic properties, except in the vinicity of the 
superfiuid transition point. Second, using the measured universal functions as a 
benchmark, we examine the applicability of a quantum cluster (virial) expansion method 
[5^ [55] . We find that for a trapped gas, up to the leading interaction effect (second 
order), the expansion is quantitatively reliable down to T ^ O.TTp. This limit can be 
decreased further to T ~ 0.4Tp, with inclusion of higher order virial coefficients (i.e., 
up to fourth order). Thus, we demonstrate clearly the usefulness of quantum cluster 
expansion in the study of a normal, but strongly interacting quantum gas. 

Finally, we note that the temperature of a trapped Fermi gas at unitarity is often 
characterized by a non-interacting temperature Tj obtained by an adiabatic sweep to the 
ideal gas limit. Using the accurate experimental data for S(T), we calculate the relation 
T{Ti), which shows an apparent kink at low temperatures. We therefore determine a 
characteristic temperature (T/Tp)o = 0.19±0.02 or (Ti/Tp)Q = 0.16±0.02 for a trapped 
Fermi gas at unitarity, below which the thermodynamic functions start to deviate from 
normal Fermi-liquid-like behavior due to pairing effects. 
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This paper is organized as follows: in Sec. 2 we briefly review different strong- 
coupling perturbation theories and the high temperature quantum virial expansion 
theory. A comparison for the bulk universal function h{C,) of a homogeneous Fermi 
gas at unitarity is presented in Sec. 3. The validity of the quantum virial expansion in 
the uniform case is discussed. In the Sec. 4, we explain how to reconstruct the trapped 
universal thermodynamic functions from /i(C), and present a systematic comparison 
of different strong-coupling theories with the accurate experimental data, for various 
thermodynamic functions. We also examine the applicability of the quantum virial 
expansion to a trapped Fermi gas at unitarity. In Sec. 5, we calculate the actual 
temperature as a function of the non-interacting temperature at the same entropy. A 
summary and outlook is given in Sec. 6. 

2. Theoretical review 

In this section we review two types of strongly interacting Fermi theories: the strong- 
coupling perturbation theories and a controllable quantum cluster expansion theory. 
We consider a two-component Fermi gas with equal spin populations. At ultracold 
temperatures (< lOOnK), the interatomic interactions between atoms with unlike spins 
can be well- described by an s-wave scattering length a^. For the case of molecule 
formation with a broad Feshbach resonance, a simple two-species model that neglects 
the molecular field is very accurate, otherwise a full three-species model is necessary 
[56| |57] . The hamiltonian of the system can then be written as, 

^ = 5Z (^k - /i) 4^Ck^ + UY1 Ck+qtCk'-q;Ck';Ck|, (2) 
k(T kk'q 

where = ^^k^/(2m) is the fermionic kinetic energy at wave number /c, and 

1 - ^ V ^ 

U Airh^s k 2ek ^ ^ 

is the bare contact interaction renormalized in terms of the s-wave scattering length a^. 

2.1. Strong- coupling perturbation theories 

We start with a brief overview of the most commonly used strong coupling theories 
for a Fermi gas at unitarity. These are approximate many-body T-matrix theories 
[58l l59] , since no exact results are known in three dimensions. Such theories typically 
go beyond BCS theory by including an infinite set of higher order Feynman diagrams. 
The diagrams included, known as the ladder sum in the particle-particle channel, are 
still not the complete set of all possible terms in perturbation theory. However, it is 
generally accepted that a ladder sum is necessary in order to include the strong pair 
fluctuations in the strongly interacting regime. This is expected to be the leading class of 
these infinite sets of diagrams [581 EHj • However, there are differences in the procedures 
used to obtain the relevant diagrams that are included. 
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Figure 1. Diagrammatic representation of T"- matrix and pair fluctuation contribution 
to the tlicrmodynamic potential. The solid line represents the single-particle Green's 
function, while the dashed line shows the bare contact interaction. Note that, these 
diagrams are valid at all temperatures, both below and above the normal-superfluid 
transition. However, below the transition, the single-particle Green's function should 
be a 2 X 2 matrix. 



In more detail, we show in Fig. la the diagrammatic structure of the T-matrix 
[SHI [59], t{Q), where the sucessive two-particle scattering between fermions with unlike 
spins is taken into account to infinite order. This forms a ladder structure, with the 
solid line and dashed line representing, respectively, the single-particle Green function 
G and the interaction U . Consequently, as an effective interaction the T-matrix can be 
diagrammatically represented by, 

t{Q) = U + UGGU + UGGUGGU + ■ ■ • , (4) 

by summing all the successive scattering process. In the normal state with contact 
interactions, the ladder sum can be conveniently calculated as, 

(5) 



(k, iWm), and q and k are wave vectors, while 



Here and throughout, Q = (q, K 
Vn = 2n7TkBT and Um = (2?t- + l)7rkBT {n = 0, ±1, ±2, ■ ■ ■) are bosonic and fermionic 
Matsubara frequencies, respectively. 

Different T-matrix theories differ in their choice of the particle-particle propagator 

xiQ), 

X (Q) = (^) (Q-K), (6) 

and the associated self-energy, 

^{K) = j:QtiQ)GAQ-K), (7) 

where we have introduced an energy-momentum summation, J2k = f^BTJ^uj,^ Sk- The 
subscripts a, /3, and 7 in the above equations may either be set to "0", indicating a 
non-interacting Green's function 

1 



Go{K) 



iUm — h ]i? /2m + /i 



(8) 
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or be absent, indicating a fully dressed interacting Green's function. In the latter case 
the Dyson equation, 



potential /i, is fixed by the number equation, = 2 lim^^o+ G(i^)e*'^™'^ [60]. By 
taking different combinations of a, /3 and 7, there are six distinct choices of the T- 
matrix, for which a notation of {GaGi3)G^ will be used. As noted earlier, there is no 
known a priori theoretical justification for which T-matrix approximation is the most 
appropriate. 

It is important to note that, while having the same diagrammatic structure, the 
T-matrix above and below the superfluid transition temperature are different, due 
to the use of different Green's functions Gq or G. In the superfluid phase below T^, 
the Green's function has to be a 2 x 2 matrix, accounting for f/(l) symmetry breaking. 
Accordingly, an additional parameter, the order parameter, appears. 

The simplest choice, {GqGq)Gq, was pioneered by Nozires and Schmitt-Rink for 
a normal interacting Fermi gas [21], with a truncated Dyson equation for the Green's 
function, i.e., 



This was shown to be equivalent to including the Gaussian pair fluctuations in the grand 
thermodynamic potential [221 El], which is shown diagrammatically in Fig. lb. The 
NSR theory was extended recently to the broken-symmetry superfluid phase by several 
authors [251 ES, [291 IHH [62]. However, some of these approaches involved additional 
assumptions to reduce computational difficulties. 

A full extension of the original idea of Nozires and Schmitt-Rink to the below 
threshold regime was reported by the present authors [2^ , with the use of a mean-field 
(2x2 matrix) BCS Green's function as "Go" in fhe thermodynamic potential in Fig. lb. 
In the following we shall refer to this extension as a Gaussian pair fluctuation or GPF 
approach. Numerical calculations were then performed at the BEC-BCS crossover for 
the equation of state of a homogeneous Fermi gas. Compared to the zero-temperature 
QMC simulation for ground state energy [12] , we found that the extended GPF approach 
works extremely well in the superfluid phase. It provides a quantitatively reliable 
description of the low-temperature thermodjTiamics of a strongly interacting fermionic 
superfluid. 

In greater detail, in our theory below the critical temperature, the contribution of 
T-matrix pair fluctuations to the thermodynamical potential takes the form (see, for 
example. Fig. lb). 




(9) 

, the chemical 




G (K) = Go (K) + Go (K) S (K) Go (K) . 



(10) 



= - V In det 




Xii (Q) X12 (Q) 

xuiQ) xii{-Q) 



(11) 
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where 




1 



Xl2 



(12) 



K 



are respectively the diagonal and off-diagonal parts of the pair propagator. Here, Qu 
and Qi2, used in Fig. lb for the single-particle line, are BCS Green's functions with a 
variational order parameter A. Together with the mean-field contribution 



f{x) = 1/(1 -|- e^/'^^"^), we obtain the full thermodynamic potential f2 = fig + 

All the thermodynamic functions, including the total energy E and total entropy 
S, can then be calculated straightforwardly following thermodynamic relations. For 
consistency, in our formalism we determine the order parameter using the gap equation 
dflo/dA = 0. Together with the number equation, n = —dQ/dn, we solve iteratively 
the two parameters /i and A. This theory with bare BCS Green functions in the 
pair propagators constitutes the simplest universal description of strongly interacting 
fermions, including the essential contribution from the low-lying collective Bogoliubov- 
Anderson modes. As such, this type of theory may have useful applications to other 
types of strongly interacting fermionic superfiuids. 

The GPF and NSR approximation does not attempt to be self-consistent. More 
sophisticated strong-coupling theories can be obtained by using dressed Green functions 
as pair propagators. For example, one may consider a {GG)G approximation, with a 
fully self-consistent propagator. This was investigated in detail by Haussmann et al. 
[23| [32] . both above and below the superfluid transition temperature. One advantage 
of the self-consistent {GG)G approximation is that the theory satisfies the so-called 
$— derivable approach to the many-body problem due to Luttinger and Ward, in which 
the exact one-particle Green functions play the role of an infinite set of variational 
parameters. The {GG)G theory is thus conserving. 

An intermediate scheme with an asymmetric form for the particle-particle 
propagator, i.e., {GGq)Gq, has been discussed in a series of papers by Levin and co- 
workers [27], based on the assumption that the treatment of fluctuations should be 
consistent with the simpler BCS theory at low temperatures. Although the {GGo)Go 
theory has been explored numerically to some extent [63], a complete numerical solution 
is difficult. A simplified version was introduced based on a decomposition of the T- 
matrix t{Q) in terms of a condensate part and a pseudogap part, leading to the so-called 
"pseudogap" crossover theory [27]. In the present comparative study, we will include 
both the pseudogap modification as well as the full {GGo)Gq theory. Due to numerical 
difficulties, we shall consider the full {GGo)Go theory in the normal phase only. 

It is clear that in the GPF approximation one omits infinite diagrams that are 
responsible for the multiparticle interactions. The fully self-consistent {GG)G theory 
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and partially self-consistent {GGq)Gq theory attempt to correct for this, by modifying 
one or more single-particle Green's function in the diagrams. However, the more crucial 
interaction vertices remain unchanged. For brevity, hereafter we shall refer to the fully 
self-consistent {GG)G theory and partially self-consistent {GGq)Gq theory as GG and 
GGq theory, respectively. 

To close the subsection, we emphasize again that there are currently no general 
grounds to decide which strong-coupling theory is the most appropriate, due to the 
absence of a small controllable interaction parameter. However, as we shall see, these 
approaches do give distinct predictions which can be tested experimentally. 



2.2. Scale-invariance and universal relation at unitarity 

In the unitarity limit, due to the infinitely large scattering length, the interatomic 
distance becomes the only relevant length scale in the problem. The internal energy 
and entropy of the system therefore scale like, 

T 



E = NepfE 
S = Nksfs 



T 



(14) 
(15) 



Tf\ 

where Tp = ep/kB is the Fermi temperature, and fp and fs are two dimensionless 
universal functions. The scaling form leads to a well-known scaling identity in free 
space: 

n = -^E, (16) 

which holds as well for any ideal, non-interacting quantum gases. To show this, 
we note that at unitarity the pressure of the gas can be readily determined from 
P = — [dE/dV]j^ g. From the expression of entropy ( IT5l) . it is easy to see that holding 
the entropy invariant is equivalent to fixing the reduced temperature T/Tp. Hence, the 
only dependence of the energy on volume is through the Fermi energy, i.e., E oc V~'^^^. 
Taking the derivative with respect to the volume, one finds that P = 2E/{3V) or 
Q = —2E/3. This simple equation relates the pressure or thermodynamic potential and 
energy for a strongly interacting Fermi gas at unitarity in the same way as for its ideal, 
noninteracting counterpart, although the energy would be quite different. 

The scaling identity Eq. ( 1T6|) follows naturally from the thermodynamic relations. 
At this point, only the GPF (NSR) approach and the fully self-consistent GG theory 
satisfy it, since in both theories we can write down a well-defined thermodynamic 
potential and then derive from it other thermodynamic quantities in a consistent way. 
Other strong-coupling theory, more or less, runs into the thermodynamic inconsistencies, 
and therefore violates Eq. f lT6|) . We note that in the superfluid phase, an ad-hoc 
renormalization of the interaction strength is required in the fully self-consistent GG 
theory, in order to obtain a gapless phonon spectrum [32]. Hence, the GG theory in the 
supefluid phase does not appear to be universal without modification. In the calculations 
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with the pseudogap theory and GGq theory, we shall obtain the energy and entropy from 
the chemical potential by integrating the thermodynamic relation n = —dQ/dfi, in order 
to satisfy Eq. ( |T6i) . The detailed procedure was outlined in Ref. [38] . 

In the experimental situation with a harmonic trapping potential V{x), potential 
energy enters into the total energy. Thus, the scaling identity should be modified 
accordingly: 

n = -^E. (17) 

This is because the potential energy {V) equals to the internal energy (3/2) / (ixP (x) 
in harmonic traps and then the internal energy is a half of the total energy. The prefactor 
of 2/3 in Eq. f|T6l) is therefore reduced by a factor of 2. 

To prove the equality, we may treat the g collection of many locally 

equilibrium uniform cells {i.e., using the local density approximation). The local force 
balance arising from the pressure P (x) and the trapping potential V (x) gives rise to, 

VP (x) + n (x) VV (x) = 0. (18) 

Taking a inner product of the above equation with x and integrating over the whole 
space, we readily obtain {V) = (3/2) / dxP (x), after the use of x ■ W (x) = 2V (x) 
for a harmonic trap and the integration by part for x- VP(x). Hence, the strongly 
interacting Fermi gas at unitary obeys the same virial theorem as for an ideal quantum 
gas. 



2.3. High temperature quantum cluster expansion 

At high temperatures there is a controllable, small parameter, given by the fugacity 
z = explfi/kBT) -C 1. This is small because the chemical potential fi diverges to — oo 
at large temperatures T. In principle, all thermodynamic properties of a interacting 
Fermi gas can be cluster expanded in powers of fugacity [SU [55], even in the strongly 
interacting limit. 

The thermodynamic potential fi^^^^ of an ideal, non-interacting uniform Fermi gas 
takes the form, 

^ = -^[^ + &^V + ... + 6«.« + ...], (19) 

where A = [27rh'^ / (mkBT)]^^'^ is the thermal wavelength and b^'' = (— l)""*"-^/?!^/^ is 
the n-th virial coefficient. We use the superscript "1" to indicate the non-interacting 
systems. While at the first glance Eq. (fT9l) may be valid at 2; < 1 only, its applicability 
is actually much wider. The expansion is meaningful for arbitrary positive values of 
fugacity through an analytic continuation across the point z = 1. It is then possible to 
show that, as expected for a non-interacting Fermi gas. 



V 



(2kBT/X^) (2/v^) t^'^ In (1 + ze~') dt . (20) 
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A622' + --- + A6„2" + --- , (21) 



In the presence of interactions, the virial coefficients are modified. The 
thermodynamic potential can instead be written as, 

n - _ 2kBT 
V ~ 

where A6„ = bn — b[l\ Our key assumption in using the quantum cluster expansion 
is that the expansion of — 1)*^^^ might be applicable near to the critical temperature, 
despite the fact that the fugacity may already be much larger than unity close to the 
superfiuid transition. It is possible to test this conjecture in either BCS or BEG hmit, 
by analytically calculating the virial coefficient A6„ and hence the radius of convergence 
of the expansion. We leave this possibility in a future study. 

At unitarity, the virial coefficients are temperature independent and are known 
up to the fourth order: A62 = 1/V2, A63 = -0.35501298, and A64 ^ 0.096 ± 0.015. 
The second virial coefficient was already known 70 years ago [51]. The third coefficient 
was calculated recently by the present authors [55]. This theoretical prediction was 
confirmed experimentally in the accurate thermodynamic measurements of Nascimbne 
et al. [19]. These recent experiments were also able to determine the fourth coefficient 
empirically [19]. With these coefficients, using = the universal function hi^) at 
high temperatures may be written as, 

MC) = 1 + ^^^c-^ + ^^3r^ + ^^^c-^ (22) 

""^^^ (2/v/^)CtV2in(l + C"ie-0cit' ^ ^ 

The above discussion for a uniform gas can be easily extended to the case with a 
harmonic trap: 14rap(r) = vnuP'r'^ j2 or, in the general case of a trap with cylindrical 
symmetry, 

VtraviA = mujlp'/2 + ujy/2 (23) 

with oj = (tol^UzY^'^ . Within the local density approximation, the thermodynamic 
potential becomes position-dependent through a local chemical potential /i(r) = fio — 
Vtrapi^) or a local fugacity z{r) = e^'^^^\ Using the fact that the virial coefficients are 
constant at unitarity, the total (integrated) thermodynamic potential QtrapifJ^OyT) = 
J drQ{r) takes the from, 

2 {kBTY 



^trap 



ihoj) 



n 



3/2 



(24) 



where zq = e^^° is the fugacity at the trap center. It is easy to see that the nth virial 
coefficient in a trap, bn,trap = bn/n^^'^, is much reduced with respect to its uniform 
counterpart. 

The thermodynamic potential of an ideal trapped Fermi gas is now given by, 

^illp = - {Wf I (hcof t' In (1 + zoe-') dt . (25) 

In analogy with the uniform case, we may define a universal function htrapiCo) = 
^trap{f^o,T)/Qlllp{fio,T), where Co = e~'^'^° is the inverse fugacity at the trap center. 
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Using the expansion at high temperatures, we find exphcitly that 




2V2Afe2Co + 8A63C0/ (3V3) + Ah 



(26) 



4a r t' In (1 + Co^^e-*) dt 



This expression will be used later on in a comparison for a trapped Fermi gas at 
unitarity. 

3. Comparisons for a uniform Fermi gas at unitarity 

We consider now the comparison between theory and experiment for a uniform Fermi gas 
at unitarity. For this purpose, we calculate the universal function h{Q using different 
perturbation theories and compare the results to the experimental measurement (Fig. 
3a in Ref. [19]) This can be easily done with the known equation of state of the uniform 
unitarity gas, that is, fi(T) = epffiiT/Tp) and E{T) = NepfEiT/Tp), where and 
fp are two dimensionless functions that depends on the reduced temperature r = T /Tp 
only. Numerically, for a fixed reduced temperature r, we calculate the inverse fugacity 
C and (C) /{Nep). We then obtain. 



3.1. Perturbation theories 

The universal functions h{Q obtained in this manner are shown in Fig. 2 for 
different perturbation theories, and are compared to the experimental data (red squares) 
[19] without any adjustable parameters. We indicate the experimentally determined 
superfluid regime by thin cross lines, where C,c — 0.042. The GG theory and the 
NSR theory above were already compared with experiment by Nascimbne et al. in 
their experimental paper [19]. The comparison presented here is much more complete. 
We find that the agreement between experiment and the three T-matrix perturbation 
theories is very good for a large range of temperatures. 

In particular, the simplest Gaussian pair fiuctuation theory gives the best 
quantitative description, though it has a possibly unphysical discontinuity at the 
superfiuid phase transition. In the Gaussian pair fiuctuation theory, the universal 
function does not have values between [Cc~]gpf — 0.15 and [Cc+]Af5i?, — 0.22. This 
non-overlap region is mostly caused by the breakdown of the original NSR theory in the 
vicinity of transition [30] , which predicts a smaller chemical potential and hence a larger 
inverse fugacity C,. In addtion, the critical inverse fugacity predicted by GPF, Q — 0-2, is 
significantly larger than the experimental observation [19], Cc — 0.042. Instead, the self- 
consistent GG theory and partially self-consistent GGq theory predict more reasonable 
values for C,ci although their agreement with the experimental data of h{C,) is worse than 



c 



exp 




r 



4 hir) 

9 r5/2 Z^;- ti/2 ln(l + C-^e-*)cit' 



(27) 
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Figure 2. (Color online). Universal thermodynamic function h{() (squares), in a 
uniform gas, compared with predictions from different strong-coupling theories: the 
GPF or NSR theory (black solid line), partially consistent GGq theory (red dashed 
line), fully self-consistent GG theory (blue dotted line), pscudogap theory (cyan dash- 
dotted line), and the BCS mean- field theory (black dash-dot-dotted line). For an ideal 
Fermi gas, /i(C) = 1, as shown by a thin dashed line. The cross region indicates the 
superfluid phase. The upper x-axis shows the temperature. 



the GPF (NSR) theory. The prediction of available quantum Monte-Carlo simulations 
is also in close agreement with the experimental data |48l|49]. Clearly, it would be useful 
to have more accurate experimental data at this point, to better understand the nature 
of the phase transition. 

On the other hand, the pseudogap theory, as a simplification of the partially self- 
consistent GGq theory, deviates significantly from the experimental data. It is therefore 
not able to capture the strong fluctuations at unitarity. However, it is certainly better 
than the BCS mean-field theory, which completely ignores pairing fluctuations. 

As the tempeature decreases to zero (T — and ( ^ 0), the universal function 
/i(C — 7- 0) — 7- where ^ = 1 + (3 is the universal parameter. Different theories predict 

different universal parameters, i.e., ^bcs — 0.59, ^gpf — 0.401, and C,gg — 0.36. Thus, 
we find that, hscsiC = 0) ~ 2.17, /igpf(C = 0) ~ 3.94, and hcciC = 0) ~ 4.65. 
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3.2. Virial expansion comparisons 

Let us now focus on the high temperature regime oi C, > 1 oi z < 1. A comparison 
of experimental data to the virial expansion in Eq. ( I22l) has already been carried out 
by Nascimbne et al. [19]. This led to the confirmation of our theoretical prediction of 
the third virial coefficient A63 ~ —0.35 as well as an experimental determination of the 
fourth virial coefficient A64 ~ 0.096 ± 0.015. The third virial calculation requires an 
exact solution of a quantum three-body problem, which is known. However, the exact 
solution of the quantum four-body problem needed for the fourth coefficient, is yet to 
be theoretically obtained. 

Here, we show that this accurate experimental data can serve as a benchmark to 
determine to what extent the virial expansion is quantitatively reliable. To be concrete, 
we shall define the criterion of "quantitative" applicability greement within 10% 

relative error for the function h{Q — 1, that is, after the non-interacting background is 
removed from the universal function h{C,). For a "qualitative" applicability, we relax 
the criterion on the relative error to 50%. 

Fig. 3 compares the virial expansion predictions (up to the fourth virial coefficient) 
for h{Q — 1 with the experimental data, using artificial 50% (a) and 10% (b) relative 
errors for comparison purposes. We are then able to estimate a critical fugacity, below 
which the n^'^-order virial expansion is either qualitatively or quantitatively valid. The 
result is tabulated in Table I, where the critical fugacities have also been converted to 
critical temperatures by using the NSR equation of state, which provides an excellent 
description of the experimental data. For the 2"*^ virial expansion that accounts for the 
leading interaction effect, we determine that for a homogeneous Fermi gas at unitarity, 
it is quantitatively and qualitatively reliable above T/Tp = 1.3 and T jTp = 0.9, 
respectively. 



Order 


2^50 




(T/Tp)50 


{T/Tf),o 


Virial2 


0.7 


0.4 


0.87 


1.30 


Virial3 


1.5 


0.7 


0.52 


0.87 


Virial4 


2.2 


1.7 


0.40 


0.48 



Table 1. Qualitative (50%) or quantitatively (10%) ranges of reliability for different 
order virial expansions in free space, as indicated by the subscript. 

To close this section, we emphasize that, though the experimental data for the 
universal function h{() are very accurate, without any prior knowledge of the bulk 
equation of state we still can not determine the temperature from the discrete data. 
However, for a trapped Fermi gas at unitarity, using the bulk h{Q we can indeed 
determine all the universal thermodynamic functions such as E(T) and S(T). This 
is discussed in detail in the next section. 
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Figure 3. (Color online). Plot of universal function h{C,) — 1 as a function of fugacity, 
compared with the experimental data (empty squares). The error bars in (a) and (b) 
indicate respectively the 50% and 10% relative errors of h{C,) — 1, in accord with the 
"qualitative" and "quantitative" criterions as described in the text. 

4. Comparisons for a trapped Fermi gas at unitarity 

Let us now turn to the experimental determination of the equation of state of a Fermi 
gas in a harmonic trap at unitarity, and a comparison of this data with theory. The basic 
idea has already been outlined by Nascimbne et al. in the Supplementary Discussion 
part of their paper [19]. 

4-1. Local density approximation 

Consider, for example, the total number of atoms, N = J drn{r). Within the local 
density approximation, we may consider the trap to be isotropic trap with a trapping 
frequency tu = (ujj_ujzy^^, without loss of generality. Using n(r) = dP (r) /9/i(r) and 
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dji (r) jdr = —mu'^r, we find that, 

roo roo 

N= c/r47rr2[9P(r)/9/i(r)] = -47r/(ma;2) / rdP{r). (28) 
Jo Jo 

Noting that P = P^^^h(Q , and integrating by parts, we obtain: 

N = ^ r drP^'^ [C (r)] h [C (r)] . (29) 

Using C i''^) = e~^^^^'^ = Co exp[ma;^r^/(2/cBT)], the integration over radius r can be 
converted to an integration over the inverse fugacity. One sees that. 



4 fkr^rV /•«^ dJ\n(C/(o) 



where / (() = i^/^/n) t^^"^ In {1 + (-'^e'^) dt. Recalhng that the Fermi energy of a 
zero-temperature trapped ideal Fermi gas is Ep = ksTp = {3NY^^huj, we may rewrite 
the above equation in dimensionless form, 



T| 12 ,^rfVln(C/Co) 



The total energy of the system may be conveniently calculated by using the scaling 
relation, Eq. f|T7|) . Thus, the total energy in a trap is given by 

poo POO 

E = -Sn = 12n / drr^P{r) = Un / drr^P^^^ [( (r)] h [( (r)] . (32) 
Jo Jo 

Converting to the variable (, we find that: 



E 72 r°° rfJln(C/Co) C 
^ = ^2TfrJ dC ^ 7'°^ n-^/(C)MC). (33) 
NEp vrV^ dC Co 

The entropy follows directly from the thermodynamic relation S = {E — Q — fioN)/T. 
Using the fact that fio/Ep = —{T/Tp) In Co, we obtain straightforwardly 

i«I = 3 TiVi; + 

The coupled equations ( 13T|) . ( l33l) and ( 134|) determine, respectively, the temperature, 
energy, and entropy of a trapped Fermi gas at unitarity as a function of the inverse 
fugacity Co- In the calculations, we discretize the integral over C and take the values of 
/i (C) solely from the experimental measured data. In this way [19], we avoid the use 
of any interpolating or fitting function to the experimental data h (C). In addition, the 
statistical error of the experimental data of h (C) is reduced. The detailed procedure for 
these numerical calculations is given in Appendix A. For convenience, we shall refer 
to this trapped equation of state, re-constructed from /i(C), as the "experimental" 
measurement or data for the equations of state of a trapped Fermi gas at unitarity. 

It is readily seen that we can calculate the theoretical prediction for the trapped 
equation of state by using exactly the same local density approximation procedure, 
combined with a theoretical universal function h{Q generated from different strong- 
coupling theories for a uniform Fermi gas. We note that our numerical procedure 
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of calculating the trapped equation of state is an average procedure integrating over 
the trap and is quite insensitive to the smoothness of h[C,). Therefore, even though 
there is a discontinuity in the theoretical universal function, as in our Gaussian pair 
fluctuation theory, we obtain a much smoother trapped equation of state. As can be 
seen from Appendix A, in that case, we simply join linearly between [Cc-]gpf ^^"^ 
[Cc+]nsr remove the discontinuity of the universal function in the non-overlap region 

([Cc-] 

4-2. Trapped universal thermodynamics: E{S) 

In previous work [11], we gave experimental evidence that any strongly interacting Fermi 
gases at unitarity has universal thermodynamics. The energy and entropy relation E{S) 
measured on ^Li and ^^K atomic clouds in three different trapping potentials all fall 
precisely on a single curve. The trapped equation of state E{S) deduced from the 
experimental data of h{() by Nascimbne et al. [19] provides an independent check of 
universality, with a much improved accuracy, in a fourth different set of experimental 
conditions. This is illustrated in Fig. 4, where we plot the new measurements using 
green circles. All the four sets of experimental data follow the theoretical prediction 
given by the simplest GPF approximation. In particular, the difference between our 
theory and the new measurement is nearly indistinguishable, as shown clearly in Fig. 
4b for the low temperature regime. This gives so far the strongest evidence for fermionic 
universality. 

To better visualize the difference between theory and experiment, we follow the 
strategy used in our previous comparative study [3H] and calculate the interaction energy 
Eint = E — EjG, which is the difference of energies between an interacting Fermi gas (E) 
and an ideal Fermi gas {Ejg) at the same entropy. Fig. 5 shows the interaction energy 
versus entropy in a harmonic trap as predicted by different strong-coupling theories 
in comparison with the experimental data reported at Duke [H] and ENS [19]. It 
is impressive that at this much reduced scale, the new experimental data obtained by 
Nascimbne et al. at ENS (green empty circles) still appear to be very smooth, suggesting 
an absolute statistical error of about O.OlNEp or a relative error of one percent in energy, 
although there may be systematic errors at this level. 

This error bar is already much smaller than the difference among different T- 
matrix approximations of GPF, GGq and GG, which is roughly of the order 0.05NEf- 
It is clear that below threshold the GPF approach provides the closest prediction to 
the new accurate measurement below threshold, with a difference at most O.OlNEp. 
However, well above threshold the fully self-consistent GG theory gives better agreement. 
We note here that the above threshold GG theory is also universal, without ad-hoc 
renormalizat ions . 
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Figure 4. (Color online). E{S) of a strongly interacting atomic Fermi gas of either ^Li 
or ''°K atoms in different trapping potentials scale into a single theoretical curve, which 
is predicted by the GPF and NSR theory. Lower panel highlights the low temperature 
regime. The cross region indicates the superfluid phase below an experimental critical 
entropy {S / NkB)c — 1-56. The determination of the critical entropy is described in 
Sec. V. The upper a::-axis in (a) plots the temperature. 



4-3. Trapped universal function hfrapiCo) 

In analogy with the comparison between theory and experiment for the universal 
function h{() for a homogeneous Fermi gas at unitarity, we consider now the comparison 
for the trapped universal function htrapiCo) = ^trap/^tlap- Because of the proportionality 
between energy and thermodynamic potential at unitarity (i.e., the scaling relation), it 
is convenient to calculate htrapiCo) using the expression, 

''"■'"^°'-9(r/r,)'/~i^in(i + cV.)di' 
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Figure 5. (Color online). Interaction energy as a function of entropy for a strongly 
interacting Fermi gas in a harmonic trap. The experimental data (symbols) are 
compared with the predictions of three different strong-coupling theories. The cross 
region indicates the superfluid phase below the critical entropy {S/NkB)c — 1.56. The 
temperature at different entropy is plotted on the upper a;-axis. 

where the denominator is simply E^^^ / [N Ep) at given inverse fugacity Cq and 
temperature T. 

Fig. 6 compares the theoretical predictions for the trapped universal function from 
different strong-coupling theories, compared with the experimental measurement at ENS 
(empty squares) [19]. Both the experimental data and the theoretical GPF (NSR) 
predictions for the trapped universal function htrapiCo) ^low much smoother after 
integrating over the trap. Thus, the discontinuity of the normal-superfluid transition in 
the GPF prediction of h{() disappears completely in htrapiCo) , due to trap-averaging. 
Comparing the different T-matrix fluctuation theories in a trap, one sees that the fully 
self-consistent GG theory gives a slightly better agreement with the experimental results 
than the GPF and GGq theories at high temperatures, with an inverse fugacity (q > 2. 
At low temperatures where Co ^ 1) however, the simplest GPF approach provides the 
best description. 

We note that, in harmonic traps the trapped universal function htrapiCo = 0) = 
^-3/2 jg exactly the same as the uniform universal function at zero temperature. 
Therefore, htrap,BcsiCo = 0) ^ 2.17, htrap,GPFiCo = 0) ~ 3.94, and htrap,GGiCo = 0) ~ 
4.65. However, these limiting values seem to be difficult to reach, compared to the 
uniform case. 
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Figure 6. (Color online). Universal function htrapiCo) in a harmonic trap. Here, Co 
is the inverse fugacity at the center of the trap. Different theoretical predictions (lines 
as indicated) are compared with the experimental measurement (squares). The cross 
region indicates the superfluid phase below a critical inverse fugacity {Co)c — 0.042. 
The upper x-axis shows the temperature. 

4-4- Quantum virial comparisons: h{() 

We now examine the applicability of the quantum virial expansion for a trapped Fermi 
gas at unitarity, by using Eq. (126|) for the trapped universal function. We use the same 
idea as in Fig. 3 and the same criterion for "qualitative" and "quantitative" reliability. 
In Fig. 7, we report the successive virial expansion as a function of fugacity up to the 4*^* 
order, compared with the experimental data. The estimates of the critical fugacity and 
of the critical temperature for different orders are tabulated in Table 11. To the leading 
second order, we find that the expansion is quantitatively reliable for temperatures down 
to T ~ 0.7Tp, much smaller than we found for a homogeneous Fermi gas at unitarity. 

This much wider applicability is due to the significantly reduced higher order 
virial coefficients in a harmonic trap, i.e., &2,tmp = ^2/(2a/2). It is readily seen that, 
with inclusion of higher order virial coefficients, the accuracy of virial expansion can 
be improved. Up to the known fourth virial coefficient, we find that the bound for 
quantitative applicability decreases further to T ^ OATp, which is a typical experimental 
temperature for a Fermi gas in its normal state. We thus show that the quantum virial 
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Figure 7. (Color online). Examination of the reliability of the quantum virial 
expansion for a trapped Fermi gas at unitarity. This is basically the same as what 
shown in Fig. 3, but now the plotted curves are computed for a trapped Fermi gas. 



Order 


^50 






(T/Tp)io 


Virial2 


1.5 


0.5 


0.45 


0.66 


VirialS 


2.7 


1.4 


0.36 


0.46 


Virial4 


3.7 


2.1 


0.33 


0.39 



Table 2. Qualitative (50%) or quantitatively (10%) ranges of reliability for different 
order virial expansions in trapping potentials, as indicated by the subscript. 



expansion method is a very useful tool for understanding the properties of a normal, 
strongly interacting Fermi gas in a harmonic trap. 
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4-5. Thermodynamic functions E{T) and S(T) 

As we mentioned earlier, in addition to the energy-entropy relation E{S), the 
measurement by Nascimbne et al. [19] was able to re-construct a complete set of 
thermodynamic functions in harmonic traps, such as E{T) and S(T). This provides us 
with a unique opportunity for a systematic comparison between theory and experiment 
for a trapped Fermi gas at unitarity, without the use of any fitting functions or adjustable 
parameters. 
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F 

Figure 8. (Color online). Universal thermodynamic function E{T). Different 
theoretical preditions (as indicated) are compared with the experimental data (empty 
squares). The cross region indicates the superfluid phase below a criticla temperature 
{T/Tp)c — 0.19. This experimental critical temperature of a trapped Fermi gas at 
unitarity is determined in Sec. V. 

Fig. 8 presents the comparison for the total energy as a function of temperature. As 
anticipated, the simplest GPF approach provides the best quantitative agreement with 
the experimental data. However, the GPF theory predicts a larger normal-superfluid 
transtion temperature, {T/Tf)c — 0.27, as indicated by a small bump. In contrast, 
the less accurate self-consistent GG theory and partially self-consistent GGq theory 
predict that {T/Tf)c — 0.21, which is much closer to the experimental observation of 
(T/T^)e^ 0.19. 

It is interesting to note that the GPF curve in the figure was first calculated by the 
present authors [30], [M] and was compared to the heat capacity measurement reported 
by Kinast et al. [13]. However, at that time, the temperature was not independently 
calibrated, due to the absence of a reliable thermometry in the strongly interacting 
regime. An empirical temperature was used, obtained by fitting the integrated one- 
dimensional density profile to an ideal Thomas-Fermi distribution. 

In this earlier comparison, empirical temperatures were converted to actual 
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temperatures using the pseudogap theory [13]. As a consequence, the resulting 
experimental data of E(T) appeared to agree well with the pseudogap theory [T3] . 
This is in sharp contrast to what is shown in Fig. 8, where the pseudogap theory 
clearly fails to account for the strong pairing fluctuations at either low temperatures 
(T < O.lTp) or high temperatures (T > O.STp). We therefore conclude that while the 
empirical temperature approach provides a rough thermometry, its model-dependence 
and insensitivity to the actual temperature makes it less useful as a tool for accurate 
comparison of theory with experiment. 




F 

Figure 9. (Color online). Universal thermodynamic functionS'(T). Different 
theoretical preditions (as indicated) are compared with the experimental data (empty 
squares). The cross region indicates the superfluid phase. 

Fig. 9 shows the comparison between theory and experiment for the entropy as 
a function of temperature. We find again that the GPF approach gives an overall 
best fit to the experimental data. Compared to the case of total energy, however, 
the effect of pairing fluctuations on the entropy is less significant. As a result, all the 
perturbation theories predict a similar entropy curve. Their difference to the ideal Fermi 
gas prediction (thin dashed line) is also small. This provides a justification for a recent 
calibration strategy used for determining the entropy of a weakly interacting Fermi gas 
[11], in which the entropy of a ^Li cloud at a magnetic field B = 1200G is assumed to 
be close to that of an ideal Fermi gas. 

4-6. Quantum virial expansion for E{T) and S{T) 

Finally, we may check the applicability of quantum virial expansion by using the 

thermodynamic functions E{T) and S{T). This is illustrated in Fig. 10, where the 

theoretical predictions of a virial expansion up to 4*^^ order for E{T) and S{T) are 
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Figure 10. (Color online). E{T) and S{T) for a trapped Fermi gas compared to the 
quantum virial expansion predictions. The cross region indicates the superfluid phase. 

compared with the experimental data. We may determine directly from the figures the 
critical temperature related to the reliability of virial expansion. These are in agreement 
with the values listed in Table 11. 

5. Thermometry of a trapped Fermi gas at unitarity 

We have noted that the temperature of a strongly interacting Fermi gas is difficult 
to measure experimentally in ultra-cold atom experiments. Unlike the situation with 
cryogenic experiments in the past, ultra-cold atoms are completely insulated by a 
high vacuum from any external reservoir at a known temperature. A useful way to 
quantify the temperature is to measure a non-interacting temperature Tj of an ideal, 
non-interacting Fermi gas. This can be easily measured from the density profile, before a 
slow, adiabatic sweep to the Feshbach resonance. Since the entropy of a non-interacting 
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Fermi gas is known, and is unchanged in an adiabatic sweep, this is essentially an entropy 
measurement. This procedure was first adopted by Regal et al. at JILA [T^. 
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Figure 11. (Color online). T is shown as a function of a non- interacting temperature 
Ti, which can be measured before the adiabatic sweep to the Feshbach resonance. 
The blue dashed line denotes equal temperatures. Inset shows clearly a kink at 
low temperature regime in an enlarged scale. By linear-fitting separately the data 
around the kink, we determine a characteristic temperature {T/Tp)o ~ 0.19 ± 0.02 
or {Ti/TF)o — 0.16 ± 0.02. The cross region indicates the experimentally determined 
superfluid phase. 



The accurate determination of the universal thermodynamic function S{T) in the 
last section presents a model-independent way to re-calibrate the Tj thermometry for a 
trapped Fermi gas at unitarity. By equating S(T) and SjciTi), where Sjg is the entropy 
of an ideal Fermi gas, we can express the temperature T of strongly interacting Fermi 
gases as a function of the non- interacting temperature Tj. The result is reported in 
Fig. 11. The inset emphasizes the normal-superfiuid transition regime. Note that the 
isentropic conversion to resonance tends to decrease the temperature so that Tj is always 
somewhat below the temperature T at unitarity. This can be seen from the dashed line 
for which T = Tj. 

We may identify three different temperature regimes from the figure. At 
temperatures Tj > 0.3Tf, the calibration curve is nearly parallel with the equal 
temperature line. To a good approximation, we find that T ~ Tj -|- 0.04Tp. Below 
Tj = 0.3Tp, the unitarity temperature seems to decrease slightly faster with decreasing 
non-interacting temperature. However, at a characteristic temperature (Tj/Tp)c — 
0.16 ± 0.02, this trend changes suddenly and we observe that the unitarity temperature 
now tends to saturate with any further decrease of the non-interacting temperature. 
This interesting feature is clearly seen in the inset, where we linearly fit the data above 
and below (Tj/Tp)c. 
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We can identify the sudden change as the deviations of thermodynamic properties 
away from normal Landau-Fermi-hquid behavior [191 HH] and therefore determine a 
characteristic temperature (T/Tp)o ~ 0.19 ± 0.02. This value agrees very well with the 
one that deduced from the homogeneous critical temperature by Nascimbne et al. [19] 
and the condensate fraction measurement by Horikoshi et al. [18]. The characteristic 
non-interacting temperature {Ti/Tp)o ^ 0.16 is also in very good agreement with the 
measurement at JILA [17| and the recent moment of inertia measurement at Innsbruck 
[65] . Using the experimental thermodynamic functions E{T) and S(T) in Figs. (8) and 
(9), we obtain {E/NEf)q ~ 0.68 and {S/NkB)o ~ 1.56. 

6. Conclusions and outlooks 

In conclusion, from the experimental data measuring a uniform universal function 
h{Q we have deduced a complete set of universal thermodynamic functions for a 
trapped Fermi gas at unitarity. These accurate experimental results provide a unique 
opportunity to test quantum many-body theories of strongly interacting Fermi gases. 
We have presented such a study by systematically comparing the theoretical predictions 
from typical strong-coupling theories with the experimental data. The comparison 
has no fitting functions or adjustable parameters. All the model approximations seem 
to be fluctuating around and not converging towards the accurate experimental data. 
We have found that the simple Gaussian pair fluctuation theory pioneered by Nozires 
and Schmitt-Rink [211 129] provides the best quantitative description for the universal 
thermodynamic properties of energy and entropy. Our comparison also includes a 
quantum virial expansion theory (or quantum cluster expansion theory) [SH [55]. We 
have investigated in detail the applicability of the expansion in the quantum degenerate 
regime. 

The experimental universal thermodynamic functions calculated in this work are 
extremely useful. For instance, the temperature dependence of entropy S{T) can be used 
to calibrate accurately the endpoint temperatures obtained from an adiabatic sweep of 
the magnetic field between the ideal and strongly interacting regimes. Therefore, by 
measuring an ideal Fermi gas temperature before the sweep and using the curve T(Tj) 
shown in Fig. 10, one can solve the troublesome thermometry problem for a strongly 
interacting Fermi gas. From these universal thermodynamic functions, we are also 
able to determine a characterstic temperature (T/Tp)o — 0.19 or (Tj/Tp)o — 0.16 for a 
trapped Fermi gas at unitarity, which is responsible for the deviations of thermodynamic 
properties away from normal Landau- Fermi-liquid behavior due to pairing effects. At 
these points, our analysis of the experimental universal function h{C) provides a new 
insight on the superbly precise experimental work of Nascimbne et al. [I9] . 

These thorough comparisons between theory and experiment provide a motivation 
for further developing the challenging many-body theory of a strongly interacting Fermi 
gases. It is impressive that the simplest Gaussian pair fluctuation approach gives such 
excellent agreement with the experimental data, especially in the below threshold regime 
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characterized by long-range superfluid order. Yet, it fails to predict the correct normal- 
superfluid transition temperature. More work is need to understand the reason for this, 
but at this stage, we feel that the GPF approximation serves as a good starting point 
for further theoretical work. Recalling that the GPF approximation includes only the 
two-body correlations (see, for example. Fig. la), a natural way to extend this may be 
to consider three-body or four-body correlations, in which three or four fermions interact 
with each other in the scattering process. The thermodynamic potential may be worked 
out with inclusion of all three-body or four-body scattering matrices. In this manner, we 
would recover correctly the equation of state predicted by the higher-order (i.e., S"*^ and 
^th qy^qy) quantum virial expansion theory at high temperatures. We believe that it will 
also lead to a more reasonable critical temperature and remove the spurious bend-back 
structure close to Tc as shown in the Gaussian pair fluctuation theory. 

On the other hand, the quantum virial expansion gives us another means for 
theoretical development, from a very different point of view. We have already shown the 
wide applicability of this expansion for the equation of state of a strongly interacting 
Fermi gas in harmonic traps, down to temperatures as low as ~ 0.4Tp. We have 
conjectured that it may be applicable down to the superfluid transition temperature, 
with inclusion of higher-order virial coefficients. Here, we can also develop a virial 
expansion for more crucial dynamical properties, such as the dynamic structure factor 
and single-particle spectral function, as measured recently at Swinburne [66] and at 
JILA [67]. The quantum virial expansion may therefore solve the troublesome problem 
of understanding a normal yet strongly interacting Fermi gas, although theoretical 
predictions beyond third order are not yet available. 

Finally, we note that there is huge interest in determining the detailed behavior of 
a strongly interacting Fermi gas near the normal-superfluid transition. Our comparative 
study, based on the most recent theoretical and experimental results, may provide useful 
insights for this future research. 
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